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APPROXIMATE  EXPANSION  FOR  FUNCTION 
THEORETIC  REPRESENTATION  OF  SOLUTIONS 
OF  THE  HELMHOLTZ  EQUATION 


MARK  D.  DUSTON 
ROBERT  P.  GILBERT 
DAVID  H.  WOOD 


BY  APPLYING  SEPARATION  OF  VARIABLES  TO  THE  HELMHOLTZ 
EQUATION,  WE  FIND  THAT  WE  WANT  TO  TRANSMUTE  SOLUTIONS  OF 


<p  +  k  2cp  =  Atp 


INTO  SOLUTIONS  OF 
d2 


dz2  if)  +  k  [1  +  £S(z)]0  =  ;*(//. 
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VUGRAPH  1 


We  are  primarily  interested  in  a  problem  from  underwater  acoustics.  By  a 
"function  theoretic  method"  we  mean  a  transmutation  of  the  Helmholtz  equation 
with  constant  coefficients  into  solutions  of  the  Helmholtz  equation  with  vari¬ 
able  coefficients.  The  Helmholtz  equation  is  a  partial  differential  equation 
governing  the  propagation  of  sound.  We  start  with  a  model  of  the  ocean  using 
the  depth  dependent  Helmholtz  equation,  and  apply  separation  of  variables.  The 
depth-dependent  part  of  the  separated  Helmholtz  equation  with  constant  coeffi¬ 
cients  satisfies  the  ordinary  differential  equation 

d2  0 

4>  +  4>  -  H  ■ 

dz 


This  represents  an  ocean  with  uniform  sound  speed.  The  depth  dependent  part 
of  the  Helmholtz  equation  with  variable  coefficients  satisfies  the  ordinary 
differential  equation: 

.2  y 

— r  +  k  [1  +  «s(z)]V>  -  A 0  , 

dz 


which  represents  an  ocean  where  the  sound  speed  is  a  function  of  the  depth 
Here  the  term  <s(z)  is  the  perturbation  of  the  index  of  refraction  of  the 
variable  coefficient  Helmholtz  equation.  The  small  parameter  e  is  a  measure 
of  the  size  of  the  perturbation.  ] 
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<M z)  -  <t>{ z)  + 


K(z , s)^(s) 


ds 


where  the  kernel  function  K(z,s)  must  be  determined.  We  substitute  the  trans¬ 
mutation  into  the  ordinary  differential  equation  for  the  separates  variable 
coefficient  Helmholtz  equation,  and  this  yields  a  partial  differential 
equation  for  the  kernel  function  K(z,s) 


~  +  k  «s(z)K  -  0 
3s 


Although  this  may  seem  to  be  stepping  backwards  ar  ter  going  through  the 
trouble  of  reducing  the  Helmholtz  equation  to  an  ordinary  differential 
equation  problem  we  will  demonstrate  a  method  of  approximating  the  kernel  b. 
iterative  integrals. 


t. 


2 


TD  7795 


We  now  invoke  our  boundary  conditions  for  the  Helmholtz  equation  with 
variable  coefficients.  We  ask  that  the  pressure  goes  to  zero  at  the  surface 
of  the  ocean  (a  pressure  release  condition),  which  is  V>(0)-0,  and  that  the 
boundary  beneath  the  ocean  and  its  bottom  at  the  depth  z^  be  rigid,  which  is 
the  condition  ip' (z^-O .  We  impose  a  similar  bottom  boundary  condition  on  the 
solutions  of  the  Helmholtz  equation  with  constant  coefficients.  Combining  this 
with  the  given  transmutation  we  obtain  two  additional  constraints  on  the 
kernel  function  K(z,s), 


as  K(z’zb> 


and 


2  K(z , z)  +  k2«s(z)  -  0 


ti 


These  two  conditions  and  the  partial  differential  equation  for  K(z,s)  are 
sufficient  to  uniquely  determine  the  kernel  function.  i_yistrlbut Ion/  ~ 
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One  method  of  solving  for  the  kernel  function  K(z,s)-  M(£,?7)  is  to 
transform  into  characteristic  coordinates  given  by 


2£  «  z  +  s  -  2zb  and  2rj  -  z  -  s 

We  then  get  a  partial  differential  equation  for  M(£,r?)  which  is 
2  M?r>  +k2£s(£  +  n  +  zb)M  -  0. 

O 

We  next  do  a  Born  expansion  of  M  in  the  parameter  k  .  This  yields  a  system  of 
integrals  where  the  first  coefficient  of  the  expansion  is  given  by  • 

2  M(1)(£,rj)  -  -f £s(a+zb)  da  - 
J0 

and  succesive  coefficients  in  the  expansion  are  given  by 


M(P+1) 


(£,»?) 


es(a+/3+zb)  da  d/9 
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VUGRAPH  5 


Inverting  coordinates  from  (£,??)  back  into  (z,s)  we  obtain  the  expansion 
of  the  kernel  function  K(z,s).  The  first  term  of  the  the  expansion  is 


m  r(z+s-2zb)/2  p(z-s)/2 

2  K  (z,s)  -  -  es(a+zb)  da  -  es(/8+zb)  d/3  , 

J0  J0 

and  the  succesive  terms  of  the  expansion  are  given  by 


K(P+1)(z,s) 


•(z-s)/2  p(z+s-2zb)/2 


K  r  (a,0)  es(a+/9+zb)  da  dfi  . 


We  have  now  obtained  explicit  formulas  for  the  power  series  expansion  of  the 
kernel  K(z,s)  in  the  parameter  «. 
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We  also  expand  in  power  series  in  the  parameter  e:  1)  the  solution  of  the 
depth  dependent  part  of  the  Helmholtz  equation  with  constant  coefficients. 
2)  the  solution  of  the  depth  dependent  part  of  the  separated  Helmholtz 
equation  with  variable  coefficients  and  3)  the  constant  X. 
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By  examiningthe  surface  boundary  condition  we  see  that  in  general  the 
transmutation  does  not  allow  the  surface  boundary  condition  for  the  separated 
Helmholtz  equation  with  constant  coefficients  to  be  identical  to  that  for  the 
problem  with  variable  coefficients  i.e.  <f>( 0)^0.  Therefore  we  define  a  new 
function 

6  (z)=  <t>( z)  -  <(>( 0) 

and  see  that  this  function  does  satisfy  the  same  boundary  conditions  as  0(z). 
We  then  expand  this  new  function  in  power  series  in  c. 
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ORDER  e: 

ei0)"  +  (k2  -  a'0))  e (0)  =  o 

ORDER  e: 

dm"  +  (k2  -  ^,0))  9n>  -  Am6l0)  =  (A{0)  -  k2)  <pm(0) 

ORDER  e2: 

0‘2'"  +  (k2  +  Ai0))  e[2)  -  An)6n)  -  At2)et0)  =  (At0)  -  k2)<p,2’(0)  +  a'VfO) 
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VUGRAPH  8 

Substituting  0(z)  into  the  ordinary  differential  equation  for  the  depth 
dependent  part  of  the  Helmholtz  equation  with  constant  coefficients  and 
isolating  the  terms  with  corresponding  powers  in  e  we  obtain  a  sequence  of 
ordinary  differential  equations  each  with  the  same  boundary  conditions  as  the 
original  problem. 

For  the  0th  order  term  it  is  well  known  that  this  homogeneous  ordinary 
differential  equation  has  a  set  of  solutions  (normalized  sine  functions)  for 
9^'  and  corresponding  eigenvalues  for  A^  ' 

We  then  take  inner  products  of,  with  the  0C  order  equation  and 
with  the  first  order  equation  and  subtract  to  get  the  expression 


^(Dpb  fl(0)(s)  ds  _  x(1)-  f  b  (A(0)-k2)^(1)(0)  S(0)( s)  ds 

J0 


which  yields  A^b^  if  we  know  the  value  of  the  constant  ^^(O) 
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We  obtain  the  value  of  the  constant  by  using  the  power  series  ex 

pansion  in  the  transmutation  aquation  and  evaluating  at  z-0 ,  this  gives 


0(1)(O) 


■i: 


bK(1)(0,s)  0(O)(s)  ds 


We  next  give  an  explicit  formula  for  A*^  in  terms  of  the  known  expressions 
for  ^(0)  and  the  expansion  of  the  kernel  function. 


(1) 


(2n-  1)tt 
-  2zb 


K  *nU<0> 


Integrating  by  parts  and  making  a  change  of  variables  we  can  show  that  the 
result  reduces  to 


\ 


(1) 


2k!  rZb 

zb  ^  0 


ts(i 


)  ^sin 


( 2m-l)7r^ 


i2 


2z, 


df 


which  is  the  same  result  given  by  Titchmarsh  in  his  classical  application  oi 
perturbation  theory. 
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Going  back  to  the  ordinary  differential  equation  for  9  and  the  first 
order  in  e  expansion,  we  must  solve  an  inhomogeneous  ordinary  differential 
equation 

0(1)„  +  (k2  _  A(0)^(l)  _  Ad),(0)  +  (A(0)  .  k2)<^(1)(0) 

for  We  solve  this  by  variation  of  parameters  and  obtain  an  explicit 
formula  for  the  coefficient  of  the  power  series  expansion.  For  general  solu¬ 
tions  to  the  homogeneous  equation  we  use 


sin- 


(2n-  1)ttz 

2zk 


and 


y2  - 


2  (2n-l)rrz 

-  cos - ^ - 

Zu  2zk 


The  solution  is  a  linear  combination  of  these  functions  and  the  particular 
solution  given  by 


il>- 1 


Zb  fl(0) (z-s) 


(A(1V0)(s)  +  (A(0)-  k2)0(1)(O))  ds 


2  (0) 

k  -  J 


is 


,(D 


The  form  of  the  solution  therefore 
boundary  condition  at  the  surface  0^(0)  __ 

mutation  preserves  the  boundary  condition  9^  '' (z^)= 
from  this  boundary  condition. 

10 


9' 

^ives 
'  ~b; 


ciy: 


+  c2y2  +  ^p1}  -  Thi- 

Because  tne  trans- 
cannot  be  determined. 
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We  substitute  the  function  0^^  into  the  equation  that  determines  0^^ 
and  use  this  with  the  transmutation  to  determine  the  first  order  in  e  coeffi¬ 
cient  of  the  power  series  expansion  of  0.  We  ask  that  this  solution  of  0  to 
first  order  in  e  be  L2  normalized.  It  is  this  normalization  constraint  that 
allows  us  to  give  a  unique,  explicit  formula  for  the  first  order  coefficient. 
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We  compare  this  transmutation  result  to  the  classical  result  or 
Titchmarsh  where  this  coefficient  is  given  only  in  terms  of  an  infinitc- 
Fourier  series. 

This  points  to  an  important  difference  between  the  two  approaches.  In 
practical  application  both  approaches  may  suffer  inaccuracy  due  to  truncation 
of  the  power  series  in  however  the  transmutation  method  does  not  suffer 
from  any  error  due  to  truncation  of  an  infinite  Fourier  series  expansion.  The 
transmutation  term  is  in  fact  the  evaluation  of  the  infinite  sum. 
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